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Abstract 

Nitrogen is the dominant impurity in the majority of natural and synthetic diamonds, and 
the family of nitrogen vacancy-type (N n V) defects are crucial in our understanding of defect 
dynamics in these diamonds. A significant gap is the lack of positive identification of N2V - , 
the dominant charge state of N2V in diamond that contains a significant concentration of 
electron donors. In this paper we employ isotopically-enriched diamond to identify the 
EPR spectrum associated with 15 N2V”and use the derived spin Hamiltonian parameters 
to identify 14 N2V - in a natural isotopic abundance sample. The electronic wavefunction 
of the N 2 V" ground state and previous lack of identification is discussed. The N2V - 
EPR spectrum intensity is shown to correlate with H 2 optical absorption over an order of 
magnitude in concentration. 

Keywords: Diamond, Nitrogen, Nitrogen Vacancy, N 2 VV H2 


1 


I. INTRODUCTION 


Nitrogen is the most common impurity found in most natural and synthetic dia¬ 
monds (with the possible exception of hydrogen), and nitrogen-containing complexes 
have attracted significant research interest during the past 50 years. A significant 
proportion of recent research into diamond has concentrated primarily on the quan¬ 
tum optoelectronic properties of the negatively-charged nitrogen-vacancy (NV _ , a 
substitutional nitrogen atom nearest-neighbor to a vacant lattice site) defect and its 
potential use in quantum metrology and sensing.-"^ NV° has been well characterized 
by optical^ and electron paramagnetic resonance (EPR) 7 techniques: the relative 
concentrations of NV° and NV~ are determined by the availability of a suitable 
electron donor (e.g. single substitutional nitrogen N s ). 

NV°/~ belong to the N n V family of defects, where n = 1-4: the defects N 2 V°, 
N 3 V° and N 4 V° are identified with the H3 (zero phonon line (ZPL) 2.465eV),-^ N3 
(2.985 eV}^ and B-center (infrared absorption band) optical bands that have been 
extensively studied.—~- 14 A notable void in the understanding of N n V-type complexes 
is the lack of identification of N 2 V~ (see figure flail : a number of candidate optical 
and paramagnetic signatures have been postulated,— 1 ^ the most promising of which 
is the optical H2 band, however no conclusive identification has been made of the 
defect by EPR. 

The H2 optical band was first reported in 1956 as a broad absorption feature ob¬ 
served after irradiation and annealing of natural diamond.— At 80 K, H2 is observed 
in both luminescence and absorption with a ZPL at 1.257eV and an accompanying 
vibronic band. Uniaxial stress measurements of the H2 ZPL assigned a symmetry 
of ,— while photochromism charge balance studies suggest that the defect is the 
negative charge state of the H3 optical defect.— 

The assignment of H2 and H3 to different charge states of the same defect presents 
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FIG. 1. (a): the structure of the N 2 V defect in diamond. Blue (light) atoms are nitrogen, 
and nominal sp -type orbitals are drawn for illustration. The defect has *^ 2 V symmetry; the 
(110) mirror plane is shown, (b): the effective one-electron picture produced by starting 
with the structure of the vacancy, lowering the symmetry and adding two (three) electrons 
for H3 (H2) using the procedure described by Lowther. 17 Electronic occupancy of orbitals 
indicated by vertical bars, with arrows indicating electronic absorption transitions between 
ground and excited states of symmetry Ai or Bi. Adapted from Lawson. 15 

a potential problem: the ground state of N 2 V~ should be EPR-active (with spin S = 
|), but has not been observed.— One candidate EPR spectrum has been tentatively 
attributed to N 2 V~, but the given spin Hamiltonian parameters are a poor fit to the 
limited experimental data.— 

The electronic structure of N 2 V°/N 2 V _ (see figure ITbl) precludes any simple spin- 
polarization mechanisms such as those seen in NV - and SiV 0 ;- however the known 
one-electron ionization behaviour— >22 from N 2 V - (paramagnetic) to N 2 V° (diamag¬ 
netic) makes the N 2 V system a candidate for long-lived nuclear spin memory pro¬ 
tocols such as those recently exploited in silicon.— 

In this paper, an EPR spectrum is identified with 15 N 2 V~ in a treated 15 N-enriched 
diamond; 15 N and 13 C spin Hamiltonian parameters are extracted from the EPR 
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spectrum and used to aid identification of 14 N 2 V~ in a sample with C & N isotopes 
in natural abundance. The spin Hamiltonian parameters are discussed and are em¬ 
ployed to illustrate the difficulty in identifying the spectrum in natural abundance 
samples. 


II. SYNTHESIS OF ISOTOPICALLY ENRICHED DIAMOND 


Isotopic enrichment is a well-established technique in the study of defects in dia¬ 
mond: nitrogen enrichment has been employed for over 30 years;— carbon enrichment 
for at least 20.— Recent reports have concentrated on carbon isotopic enrichment in 
chemical vapor deposition (CVD)-grown diamond; however CVD synthesis becomes 
problematic with high levels of gas-phase N 2 and hence high pressure high temper¬ 
ature (HPHT) is still preferred when doping with ^ 10 ppm nitrogen in the solid 
phase. 

Atmospheric nitrogen adsorbed (absorbed) onto (into) the growth source and cap¬ 
sule materials is readily incorporated into the synthesized diamond ([N s °] ^ 100 ppm 
are typical) if preventative measures are not employed. Such measures include chem¬ 
ical nitrogen traps (so-called “getters”) added to the growth materials, or outgassing 
all growth materials pre-synthesis and sealing the growth capsule. The sample used 
in this experiment was grown by evacuating and outgassing a sealed HPHT syn¬ 
thesis capsule and subsequently backfilling with 15 N-enriched N 2 gas:— ^ 95 % of 
incorporated nitrogen was 15 N — nominally identical to the source gas. 
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III. EXPERIMENTAL DETAIL 


A. Sample 

Both N s ° and A-centers ((N s — N s )°) are effective traps for mobile vacancies. 
If N s ° is the most abundant impurity post-irradiation (to produce vacancies and 
interstitials; we will not concern ourselves here with the interaction of self-interstitials 
with nitrogen defects nor the charge state of the vacancy, as neither has a bearing on 
the following discussion), then upon annealing to temperatures where the vacancy 
is mobile 600 °C) the dominant aggregation mechanisms will be N s ° + V —> 
NV° and N s ° + NV° —> N s + + NV - ; similarly if A-centers are abundant then the 
dominant process becomes (N s — N s )° + V —y N 2 V 0 . If both A-centers and N s ° are 
present they will both trap mobile vacancies and N s ° will donate charge to N 2 V, 
undergoing the additional process N s ° + N 2 V° —» N s + + N 2 V - . NV is stable to 
approximately 1500 °C; at this temperature the reaction N s ° + NV° —>• N 2 V° and of 
course N S °+N 2 V° —> N S + +N 2 V - start to become significant. N 2 V is itself only stable 
to approximately 1600 °C and anneals out with the production of A-centers (N 2 V —> 
(N s — N s ) T V), where the vacancy is recycled to promote further nitrogen aggregation 
/ migration.— Thus when irradiated type lb (N s ° is the dominant impurity in starting 
material) diamond is annealed, only a very narrow window exists where N 2 V is 
produced over NV — it is challenging to produce large concentrations of N 2 V by 
this route. However, if the type lb diamond is subjected to HPHT treatment to 
convert a substantial fraction of N s ° (but not all) to A-centers before irradiation and 
subsequent annealing at 800 °C, then significant concentrations of N 2 V - (and N 2 V 0 ) 
can be produced. 

The sample used in this study was annealed using the above protocol: it first 
underwent high temperature (^ 1900 °C) and pressure (of order 5 GPa) treatment 


5 


for 1 h to produce nearest-neighbor nitrogen aggregates. The sample was then elec¬ 
tron irradiated with 4.5 MeV electrons to create vacancies. This processing pro¬ 
duced approximately [V°] = 5(1) ppm, and [( 15 N S — 15 N s )°] = 45(3) ppm, [ 15 N°] = 
10(1) ppm and [ 15 N+] = 6(2) ppm, accounting for pre-treatment concentration of 
[ 15 N°] = 105(5) ppm. Finally, the sample was further annealed at 800 °C for 14 h to 
produce nitrogen vacancy aggregates. Final paramagnetic defect concentrations were 
approximately [N s °] = 5.0(4) ppm, [NV“] = 1.6(3) ppm and [N 2 V~] = 1.8(2) ppm: 
the sample was notably photochromic and hence concentrations were subject to re¬ 
cent sample history. An NIR-visible absorption spectrum of the sample is given in 
figure [2al 




Energy (eV) 


N 2 V‘ Concentration from EPR (xl0 17 cm' 3 ) 


(a) 


(b) 


FIG. 2. (a): NIR-visible absorption spectrum of the sample following high temperature 
high pressure annealing, electron irradiation and further annealing. Strong absorption 
bands include H2 (1.257eV) and NV _ (1.945eV). (b): correlation between the integrated 
absorption of the H2 optical band and the 15 N 2 V _ EPR intensity in neutron-irradiated and 
annealed samples. 
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IV. RESULTS 


A ^ 2 v defect such as N 2 V may be aligned along one of six equivalent orientations, 
each defined by its principal {110} plane (see Supplemental Material for illustra¬ 
tions of different orientations; all following directions are for the orientation given in 
figure [Tal) . We expect the g-tensor to be anisotropic with principal directions along 
[110], [110] and [0 01]. When observed in EPR with the Zeeman held B along a 
(0 01) direction, the (^-anisotropy splits the orientations into two distinct groups: two 
orientations with their defining (110) axis perpendicular to B ; and four orientations 
with their axis at 45°. N 2 V - contains two equivalent nitrogen nuclei - in this case 
both nuclei are 15 N and therefore possess I — |. A reasonable assumption is that 
the 15 N hyperhne interaction is approximately along the (111) direction connect¬ 
ing each nitrogen to the vacancy, and hence all nitrogen hyperhne interactions are 
equivalent when R||(0 01). The expected structure for an EPR spectrum recorded 
with 51| (0 01) is therefore two sets of lines (with intensity 2:1 according to the 4:2 
orientation ratio) each split by two equivalent I = \ nuclei. 

The experimental B\\ (0 0 1) spectrum given in hgure [3] displays the expected struc¬ 
ture: two sets of lines with relative intensities 1:2:1. Similar arguments apply for both 
B || (111) and R||(l 1 0) spectra. These spectra unambiguously identify an S = | de¬ 
fect containing two equivalent I — \ nuclei at near 100 % abundance, indicating 
15 N. 

NIR-visible optical absorption measurements were performed at room tempera¬ 
ture on several diamond samples treated with neutron irradiation and subsequent 
annealing. In all samples the intensity of the N 2 V _ EPR spectrum was found to 
correlate with the integrated absorption of the H2 ZPL — see hgure [2b] This result 
gives additional evidence to the assignment of the H2 optical band to the N 2 V~ 
defect.— 1 ^ 
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FIG. 3. A comparison of experimental (a) and calculated spectra (b) at a microwave 
frequency of approximately 9.755 GHz. Calculated spectra generated with EasySpin2£ using 
the spin Hamiltonian parameters determined by fitting the line positions (see table H|; EPR 
linewidth was adjusted to fit the experimental data. The orientation of the magnetic field 
for each spectrum is given in the center of the figure. The calculated spectra include 15 N 2 V - 
and 14 N s ° (visible at 348.82 mT), and have been referenced to N s ° assuming that this defect 
possesses an isotropic g-value: 2.0024. 





























V. DISCUSSION 


A. The nitrogen hyperfine interaction 

In the simple case, the hyperfine interaction can be described as the sum of an 
isotropic component arising from non-zero spin density at the nucleus (Fermi con¬ 
tact), and an anisotropic dipole-dipole component. Typically the sign of the interac¬ 
tion is ambiguous: here, the sign of the 14 N and thus 15 N hyperfine was determined 
using knowledge of the sign of the quadrupole interaction (section fV Cl) . 

The isotropic component a = + A 2 + A^)/2> = +4.02 MHz (see table Q]) is 

small and opposite in sign to that expected for a negative nuclear magnetic moment. 
This implies a negative unpaired electron probability density at the nitrogen nu¬ 
cleus, and hence indicates an indirect interaction (configuration interaction, exchange 
polarization):— a similar situation was found for the ground state of NV _ , where the 
unpaired electron probability density is localized on three carbon neighbors.— For 
N 2 V" this suggests, as is confirmed by the 13 C hyperfine interaction (section IV BD . 
that the unpaired electron probability density is localized on the two carbon neigh¬ 
bors of N 2 V _ , polarizing the core states of the nitrogen and, since the nuclear mag¬ 
neton for 15 N is negative, yielding a positive Fermi contact term. Given the near¬ 
zero localization of the unpaired electron probability density on the nitrogen, the 
anisotropic dipole-dipole interaction must originate due to an interaction between 
the nitrogen and the unpaired electron probability density localized in an orbital on 
another atom. 

A crude model was employed to aid interpretation of 15 N anisotropic hyperfine 
interaction. In N 2 V°, each nitrogen atom possesses two electrons in its lone pair 
pointing into the vacancy (each nitrogen is back-bonded to three carbons; see fig¬ 
ured]), whereas both carbon atoms neighboring the vacancy possess only one electron 
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in the non-bonded orbital pointing towards the vacancy. These two orbitals interact 
to form an extended bonding orbital, which is thus fully occupied in N 2 V°. It was 
therefore assumed that the unpaired electron probability density in N 2 V~ is local¬ 
ized in an antibonding orbital formed between these two carbon atoms. An axial 
dipole-dipole interaction matrix of the form b , b, —2 b was constructed for each of the 
unpaired electron probability density locations, then each matrix was transformed 
into the crystal axes and summed to yield a macroscopic interaction matrix (the form 
measured in experiment). For each nitrogen atom the interaction was restricted to 
the {110} plane containing the nitrogen, vacancy and the appropriate next-nearest 
neighbor carbon. No geometric relaxation was included in this model, and the un¬ 
paired electron probability density was taken to be localized at two points only. 

A least-squares fit was performed against the experimental values, with the only 
free parameters being 9, the in-plane angle of the interaction and &g t , the strength of 
the interaction in MHz. 6 was measured from the (0 11) direction between the nitro¬ 
gen and the next-nearest neighbor carbon and gives the direction of the interaction; 
feat indicates the strength of the interaction (and hence dipole-dipole proximity). A 
value of 9 = 35.3° would indicate that the unpaired electron probability density is 
localized on the axis connecting each nitrogen atom to the vacancy, whereas a value 
of 9 = 0° suggests localization on the axis connecting each nitrogen atom to its 
next-nearest neighbor carbon atom. 

The best fit to the dipolar interaction is given in table [TT] and was achieved at an 
angle of 9 = 0° and b^ t = —0.46 MHz. The resulting macroscopic interaction differs 
in orientation from experiment by 3°. A simple electron- 15 N dipolar calculation 28 
yields this value for b^ t with a 15 N to carbon atom separation of approximately 
2.1 A. Reassuringly, this distance is only 18 % smaller than the next-nearest neighbor 
separation of 2.52 A in diamond. For a more precise calculation, information about 
shape and amplitude of the electron wavefunction is required in addition to knowledge 
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about the minimum energy geometrical configuration. 


B. 13 C hyperfine interaction and localization of the unpaired electron prob¬ 
ability density 

The unpaired electron wavefunction Tj may written as a summation over all atoms 
where the unpaired electron probability density is non-zero: 

N N 

Vj = ^2 71n ^ n where Vn = 1 • 

n n 

For carbon and nitrogen atoms in diamond, each atomic wavefunction consists of 
contributions from s- and p-type wavefunctions, with 

i>n = a n (p 2 s + Pnfap and a 2 n + (3^ = 1 . 

Using standard tables and methods for hyperfine interpretation,— the measured 13 C 
hyperfine values a = 240.7MHz and b = 38.4 MHz yield o? G = 15%, /3q = 85% and 
an unpaired electron probability density on each nearest-neighbor carbon atom of 
r/c = 42 %. 

The calculations above are consistent with one another and form a complete pic¬ 
ture: the unpaired electron probability density is highly localized on the two carbon 
atoms, and the nitrogen interaction is weak due to virtually zero local unpaired elec¬ 
tron probability density. Approximately 85 % of the unpaired electron probability 
density can be accounted for by the two carbon atoms nearest-neighbor to the va¬ 
cancy in N 2 V _ . This is very similar to the NV - system, where 84% of the electron 
spin density can be accounted for by the three nearest-neighbor carbons, 27 with an 
.sy>hybridization ratio of A 2 = 6.2(2) versus 5.6(2) here; these values suggest that the 
relaxation of the carbons away from the vacancy is similar in both defects. 
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C. The 14 N quadrupole interaction 


As discussed in the introduction, has not previously been identihed in 14 N- 

doped diamond. However, its identification in 15 N-doped diamond provides a route 
to identification in natural abundance samples: by scaling the hyperhne parameters 
by the ratio of the isotopic nuclear g-value,— the only unknown parameter required 
to fit a potential 14 N 2 V _ spectrum is the quadrupolar interaction with 14 N. Using 
an electron localization in the 14 N 2 p orbital of 0%, as observed in experiment, 
a quadrupolar interaction strength of Pm m —5.2 MHz can be estimated following 
previous discussions on the magnitude of 14 N quadrupolar interaction strengths in 
diamond.— Additionally, following other nitrogen- and vacancy-containing defects in 
diamond, the quadrupolar interaction can be assumed to be aligned along (111). 

An HPHT synthetic sample with natural nitrogen isotope abundance was prepared 
using the treatment procedure described in dlllL subsequent measurements recorded 
the experimental spectrum shown in figure [4l Simultaneous fitting of (0 0 1) and 
(110) spectra yielded a best-fit value of P\\ = —5.0 MHz, in close agreement with 
the estimate based on other N n V defects in diamond. The characteristic 14 N hyper¬ 
hne structure is not immediately apparent due to the large quadrupolar interaction 
leading to mixing of nuclear spin states (81 possible transitions per defect versus 16 
per defect for 15 N), with the effect that the spectrum is “smeared out”. Additionally, 
the EPR spectra of 14 N 2 V - and 14 N S ° significantly overlap (and hence obscure one 
another — see figure |4j) but are spectrally separated in 15 N-doped material (due to 
/ = 1 —y I = | for 14 N and 15 N, respectively). A combination of the inherently 
more complex spectrum and the difficulty of producing high concentrations of N 2 V _ 
without higher concentrations of N s ° explains why 14 N 2 V~ has not been previously 
identihed. 
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FIG. 4. (a): experimental B \\(0 0 1) EPR spectrum of a 14 N-doped diamond post-irradiation 
and annealing; calculated spectrum of 14 N S ° and an additional line at g = 2.002 69; dif¬ 
ference spectrum, (b): difference spectrum (as bottom left); 14 N 2 V _ spectrum generated 
using the Hamiltonian parameters given in table QT| 

VI. CONCLUSION 

An 5 = 1 spectrum containing two near-100 % abundant I — \ nuclei has been 
observed in a 15 N-doped synthetic diamond, indicating 15 N as the nuclei involved. 
The corresponding defect has symmetry. A full description of the unpaired 
electron probability density localization has been proposed and is entirely consistent 
with all expected electronic and spectral attributes of N 2 V _ , a rhombic defect in the 
N n V family of diamond defects. The vast majority of unpaired electron probability 
density in the defect is distributed over the nearest-neighbor carbon atoms of the 
vacancy, as is the case in other vacancy-type defects in diamond. 27 ' 31 
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The fitted spin Hamiltonian parameters were employed to identify the correspond¬ 
ing spectrum in an 14 N-doped synthetic diamond, where the only free spectral param¬ 
eter was the quadrupolar interaction strength — the nuclei are therefore identified 
unambiguously with nitrogen. The large quadrupole interaction and overlap of the 
14 N 2 V~ and 14 N S ° EPR spectra explains why the defect has not previously been 
identified. 

The integrated intensity of the N 2 V _ EPR spectrum has been shown to correlate 
with the integrated absorption of the H2 optical band over an order of magnitude in 
concentration, further strengthening the assignment of H2 to the N 2 V - defect. 

Isotopic enrichment has proved instrumental in the understanding of N 2 V, and 
demonstrates the importance of complementary synthesis techniques in developing a 
holistic view of the possibilities of defect engineering in diamond. Isotopic enrichment 
may be key in the exploitation of this defect as a photochromic memory with EPR 
readout.— 
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Parameter 

Value 

Direction 

9i 

2.003 45(5) [110] 

92 

2.002 74(5) [0 01] 

93 

2.002 71(5) [llO] 

Parameter 

Value (MHz) 

Direction 

15 N A 1 

+3.47(2) 

-3.5(5)° from [112] 

-A 2 

+4.51(2) 

-3.5(5)° from [111] 

A3 

+4.09(2) 

[1T0] 

13 C A„ 

+317.5(5) 

2.0(5)° from [1 11 ] 

Au 

+202.3(5) 

[110] 

14 N P,| 

-5.0(1) 

[HI] 


TABLE I. The measured spin Hamiltonian parameters for the N 2 V - defect in 15 N- and 
14 N-doped diamond. The sign of the nitrogen hyperfine interaction is determined by the 
negative quadrupolar interaction for 14 N; the sign of the carbon hyperfine must be positive 
to account for the observed electron probability density localization, ^-values referenced 
to N s ° assuming isotropic g = 2.0024. Directions relate to the defect orientation shown in 
figure [T3 — the hyperfine interaction for the second nitrogen / carbon can be generated by 
a C 2 rotation about [0 01]. Angles are given as the acute angle between the interaction axis 
and [0 01]. 
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Magnitude (MHz) 

Direction 

Calculation 

-0.58,+0.46,+0.12 

[112], [111], [llo] 

Experiment 

-0.55,+0.49,+0.07 

3.5° from [112], 3.5° from [111], [110] 


TABLE II. Results of calculations designed to minimize the difference between the calcu¬ 
lated and observed hyperfine values. The values were calculated at an angle of 9 = 0 (see 
text for details). Values can be compared to the experimental values given in table Q] by 
adding the isotropic 4.02 MHz component. 


18 





